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Abstract 

We investigate the supersymmetric electroweak corrections to the decay widths 
of the CP-odd and the heavy CP-even Higgs bosons into chargino or neutralino 
pair in the framework of the Minimal Supersymmetric Standard Model. The cor- 
rections involve the contributions of the order 0{aew'm-^i,-f/'m,'yy), 0{aew'm-t(i,^/n^w) 
and 0{af.^mi(j,'^ I m\Y) ■ The detailed calculations of the electroweak corrections to the 
following decay processes: AP / XiXi and /H^ X2X2 ^.re presented in this 
paper. We find that these relative corrections maybe rather large quantitatively, and 
can exceed 10% in some regions of parameter space. The corrections to the decay 
/H^ X1X2 <^a,n be obtained analogously, but our results show that they are very 
small and can be neglected. 
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1 Introduction 

In the Standard Model (SM) there exist quadraticahy divergent contributions to the 
corrections to the Higgs boson mass. This is the so-cahed naturalness problem of the SM. 
A good way to solve this problem is to consider supersymmetric (SUSY) extensions of 
the SM. Then the quadratic divergences can be cancelled by loop diagrams involving the 
supersymmetric partners of the SM particles exactly. The most attractive supersymmetric 
extension of the SM is the Minimal Supersymmetric Standard Model (MSSM) H, ^. In 
the MSSM, we need two Higgs doublets Hi and H2 to give masses to up and down type 
fermions and to assure the cancellation of anomalies. The Higgs spectrum consists of 
three neutral Higgs bosons: one CP-odd particle {A^), two CP-even particles {h? and 
H^), and a charged Higgs boson (H^). At the tree level, the neutral Higgs boson masses 
are determined by two input parameters (see Eq.( 2.13| )). In this paper, we choose the 



mass of the CP-odd Higgs boson rriA and tan f3 as the two input parameters in Higgs 
sector, which is defined by tan/3 = V2/V1, where vi and V2 are the vacuum expectation 
values of the two Higgs boson fields Hi and H2, respectively. The radiative corrections to 
the Higgs boson masses have been discussed in Ref.Q, ||, ^. We find that the magnitude 
of the corrections to niH and nij^ are less than 1% with our parameter choice, therefore 
they will not be considered in our calculation. 

The neutral Higgs boson decays can be divided into two types of decay modes, i.e. 
the SM decay modes and the SUSY decay modes. The first type of decay modes (e.g. 
A^/H^ T^T^, H^ 77 and ZZ or ZZ* 41) have been already carefully discussed, 
especially in the parameter space with large SUSY particle masses where the Higgs boson 



decays into sparticles are kinematically forbidden |^ . But for the second type of decay 
modes, it has been found that the decays of the MSSM Higgs bosons into charginos or 
neutrahnos may have substantial branching fractions when kinematically allowed, and 
can sometimes even dominate the usual SM decay modes |9|, |T^. Moreover, the one-loop 
calculation of the Higgs decay into the lightest neutralino pair (xiXi) is important for the 



dark matter searchjll, 12 1. Therefore, the widths of these SUSY decay modes should be 
calculated more precisely. In this paper, we will study the radiative corrections to the 
decays of and into charginos or neutralinos in detail. 

Considering the Yukawa corrections of the order 0(ae«,m^^j^^/m^), 0(ae«,m^(-fj^/rn,^) 
and 0{aew'm't{b)/^w)i we need calculate the one-loop Feynman diagrams involving the 
third generation of quarks and squarks which are the main contributions to the electroweak 
radiative corrections. The structure of this paper is as follows: in Section 2, we review 
briefly chargino-neutralino, squark and Higgs sectors of the MSSM, and give some results 
at the tree-level, in Section 3, we present our renormalization procedure, especially for 
the Higgs sector, and give some analytical results of the corrections to the decay widths, 
the numerical results and conclusions are presented in section 4, and in the Appendix, 
we present the explicit expressions of form factors for these decays and some necessary 
formulas. 

2 Related theories at tree-level 
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2.1 Chargino and Neutralino sector 

The tree-level chargino and neutralino mass matrices X, Y are 

M ^/2mw sin /? 



^ 

\/2mw cos (3 /i 



Y 



/ M' —mzsw cos (3 mzsws^''^P \ 

M mzcw cos (3 —mzcw sin (3 

—mzsw cos f3 mzcw cos (3 — 

\ mzsw sin f3 —mzcw sin j3 —fi / 



(2.1) 



To obtain the tree-level masses of chargino and neutralino, we introduce three unitary 
matrices U, V and to diagonalize X and Y: 

U*XV'^ = Md = diag(m-+,m-+), (0 < m-+ < m-+) 

Al A2 Al A2 

7\r*yAr^-'^ = Nd = diag(m^o, m^o, m^o, m^o), (0 < m^o < m^o < m^o < m^o).(2.2) 
where m-+, (i = 1,2) and m^o, (j = 1, ...,4) are the masses of charginos and neutralinos 

Ai Aj 

respectively. 

2.2 Squark sector 

The tree-level stop and sbottom mass matrices are 



M 



M? + + m| cos 2/5 (i - - cot /3) 

I + mf + |m| cos2/^a^^ 



* I mt{At- cot /?) M? + + |m| cos 2/35^ 



^ I rubiAb- fi tan /?) M? + - |m| cos 2/3s^ I 

where Mq, M^, Mj=^ are the soft-SUSY-breaking masses and Af^i, the corresponding soft- 
SUSY-breaking trilinear coupling parameters. They can be diagonalized by unitary ma- 
trices IZ^ via 

n^Mln^^ = diag(m?^ , m|) (2.4) 



where 



cos t>a sm ( 



Sm Oq COS tlq 



(-7r/2 < Bq < Tr/2, q = t, 



(2.5) 



2.3 Higgs sector 



For the convenience of the Higgs potential, we introduce two SU (2) Higgs doublets Hl, Hh 
defined as 



H2 = Hl, 



Hi = H% = iT2H 



R 



(2.6) 



where T2 is the second Pauli matrix. The Higgs potential can be divided into four parts 



(2.7) 



which represent the linear, quadratic, cube and quartic terms respectively. In terms of the 
shifted fields by [| § 



Hl 



VL/V2 + H^, + iHl, 



the linear term can be expressed as 



Hr 



vn/V2 + H^, + iHj,, 



(2i 



V, 



(1) 



H 



LO 



mLH+ I M I +- 



+ V2vrH^o 



VR 



-.2 _L ^'2 



2 I i2 9 ~^ 9 / 2 2 \ v^-i^ 

mRH+\fJ'\ + z (^iJ - ^'l) - "T- — 



(2.9) 



where miH^t^RH and m are the soft-SUSY-breaking Higgs sector mass parameters. We 
define the coefficients of and in Eq. (|2.9D to be 



Tl 
Tr 



VL 



VR 



2 ii \2 , 9 ^ 9 I 2 2\ 



2 1I 

2 ii |2, 5+5/2 2\ 2'^^ 



2'^ 
VL 

2}^ 
VR 



(2.10) 



Then the hnear term of Higgs potential can be written as 



y« = V2HSqTl + V2H^qTr. (2.11) 
The quadratic Higgs potential vjj^ in Eq. ( |2.7| ) consists of two parts which are relevant 



to charged Higgs and neutral Higgs bosons, respectively. In this paper, we only consider 
the radiative corrections to the MSSM neutral Higgs boson decays. The neutral Higgs 
part of can be written in the form Q 



y(2)ncn _ ( rrl rrl W 1^ + ^ \{ 



+ ( H,,H^o )[ _ ^^^^^ if + -'^ + ^4 j U^o J 



By defining 



^(^io) ^ ^^^^ 



Hlo \ _ crp\ ( \ - ( sin/? -cos^ \ ( 

cos/3 sin/3 y \ —A^ 



tan P = —, tan 2a = — ^— — f - tan 2/? 



vlu/j sin/3cos/3 

1 



^{m\ + m|)2 - 4m|m^ cos2(2^)^ (2.13) 



T// = sin a + Tr cos a 



Tl cos a — Tr sin a, 



(2.14) 
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we can get 



with 



hHH = [Tff (cos^ a sin /3 + sin^ a cos /3) + T/i sin a cos a sin(a — /?)] 

mw sm 2p 

bhh = : — TT^ [Th sin a cos a cos(a — /3) + T)i(cos a cos /? — sin a sin f3)] 

m\Y sin 2p 

bnh = [Tff sin(a - P) +Th cos(a - /?)] 

2mvK sm 2jJ 

boG = -^[THCOs{a-p)-ThSm{a-P)] 

Baa = -[Tfffsin^ /3cosa + cos"^ /3sina) + T/i(cos^ /3cos a — sin^ /Ssina)] 

mw sin 2p 

bcA = -^[THSm{a-(3)+ThCOs{a-P)] (2.16) 
2.4 Tree- level results 



The relevant part of the MSSM Lagrangian at the tree-level to the decays concerned in 
this paper, can be expressed as 

2 

£ = J2^-9H'^t{PLjf,* + PRj!J)xj~gh^^tiPLJ^r + PR4)xt 

+ igG^'xtiPLjf* - PRjS)xt + igA^xliPLJt* - PR4)xt] 

4 



+ ^ E [gH'^'iPLFi^* + PRP^^)x'j + gh'-^iPLF^* + PrfI;)x'j 

- igG^tiPLFS* - PrF,^)x', - igA^f,{PLF,f - PRF,j)x^] (2.17) 
where Pl,r = ^(1 T Ts); and the matrices j^'^^''^'^ ^ jp^'^^''^'^ gj-Q defined in Appendix A. 



By using above Lagrangian, we can easily obtain 



MtiH'^xtxJ 



I MM' - 1^ 



(I 4 I' + I 4 1')"^^ - (I 4"^X+ - •^jt"^;^ + 



(I -/i^ I' + I 1')"^!^ - (I Jljrn^f + 



25^ 



2/ 



(2.18) 



The tree-level decay widths can be obtained immediately after the integration over the 
phase space. 

3 Renormalization and radiative corrections 



The Feynman diagrams involving the third generation of quarks and squarks are shown 
in Fig.l (Fig.l.c and Fig.l.f). In our calculation, we use the t'Hooft gauge and adopt 



the dimensional reduction (DR) [13| scheme, which preserves supersymmetry at least 
at the one- loop level. The complete on-mass-shell scheme jl4, 15, ^] is used in doing 
renormalization. 
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The relevant renormalization constants in our calculation are defined as 

GO \ ( 1 + \bZGG ^SZga \ ( 

J ^ \ \5Zag 1 + \5Zaa )\A^ 

\ / 1 + '^6Zhh \^Zhk \(H^ 

J ^ \ \5ZhH 1 + k^Zhh )\h^ 

g ^ g + 6g, (3^p + 6(3 

U U + dU, V^V + SV, N^N + 6N. (3.1) 

All these renormalization constants can be fixed by the on-mass-shell renormalization 
conditions: 

(1) The vacuum expectation values (VEV) of and are zero ||5|: 

(0 I I 0) = (0 I /i° I 0) = (3.2) 

At the tree-level, this means Th = = 0, and therefore bHu,Hh,hh,GG,GA,AA = 0. Then 
from Eq.( 2.15| ), we can get the counterterm of V^^ 



8V^^ = 6ThH^ + dThh° (3.3) 



When we consider the radiative corrections at one- loop level, the vacuum conditions (3.2) 
tell us that 

if^ = iT" - i6TH = 



9 



where iV^'^ are renormalized one-point Green functions (renormalized tadpoles) of 
and , and iT^'^ ^—iSTjj ^ Eire the contributions from the one-loop diagrams (tadpole 
diagrams) and conterterm, respectively. 

(2) On sheU conditions for MSSM Higgs bosons [§ |]: 
Analogously, we can also get the counterterm of V^'^^^" from Eq.( |2.15 ). Furthermore, all 



the renormalized one-particle-irreducible two-point Green functions of the MSSM neutral 
Higgs bosons can be obtained directly. The expressions of them are given 

in Appendix B. In our renormalization scheme, we choose m^, mz as the input parameters 
for Higgs sector and set mA-,fnz to be the physical masses of and Z^: 

rriA = mF^^^' mz = mF^^^' (3.5) 

The relevant on-shell renormalization conditions for Higgs sector are 



■f^HHt phys.2N --(^Hhf phys.2N n (^'^) I ■ /o c\ 

zT {my )=ir {my ) = |p2=„Phys.2 = i (3.6) 

where m^^^' is the physical mass of H^. Then we obtain 
5m\ = J:^^{m\)-5bAA 



m 



P|^y-2 ^ + 5^2^ _ S^^(m2,) + 6bHH 



dZAA \p2=m\ OZhH \v^=v^l 



2 



5Zga = — [5bGA-Ty''{m\)] 
m\ 



2 



bZhu = — ^[5hHh - S^'^(m^)] (3.7) 
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with the definition 



H 



d(3 



(3.8) 



(3) For the renormaUzation of the parameters /? and a, we use the relation between 
5P and 5Zga 



6Zga = sin 2/3 



(5tan/3 Svr 5vl 



tan j3 



and impose 5vl/vl = Svr/vr, then we have 

4 



6/3 



5a = sin4Q! 



+ 



(3.9) 



sin 4/3 ' 2{m\ — m\) 
The renomahzation of the rest parameters and fields can be found in Ref . p^, p^ , p!8| , 19 1. 
By substituting Eq.(3.1) into the Lagrangian ( p. 17] ), we can obtain the counterterms 



of these vertices (see Eq.(5^) in Appendix A). Then the renormalized amplitudes can be 
expressed as 



(3.10) 



where Mc and ^Av represent the counterterm and the vertex corrections respectively. 
The divergence in A^^, can be cancelled by that in A^c exactly. It can be checked 
both analytically and numerically. The form factors of the renormalized amplitudes 
Mtot{A^/H^ xtxj) and Mtot{A'^/H'^ ^ XiX^) can be written as 



J, 



(tot)A,H 



jA,H r{v)A,H j{c)A,H 
ij ~^ ij ~'~ ij 



p(tot)A,H ^ pA,H p(v)A,H Tp{c)A,H 
ij ij "T ij "T ij 



(3.11) 



11 



where J^j^^'^ and F^J^^'^ are the terms contributed by the vertex corrections. The ex- 
pressions of these form factors are hsted in Appendix A. Then we can get the renormahzed 
decay width Ftot directly. To describe the magnitude of the radiative corrections to the 
decay width, we define the relative correction 6 = (Ftot — rtree)/rtree- This quantity will 
be used in the next section. 

4 Numerical results 

In this section, we present some numerical results for the quark-squark loop corrections to 
the chargino and nuetralino decays of and H^. Since we find the correction to the decay 
/ X1X2 very small and can be neglected, we present the results only for the 
decays of AP / xtxi /X2X2 here. In our numerical calculation, the SM parameters 
are taken to be mt = 174.3 GeV, rub = 4.3 GeV, mz = 91.188 GeV, mw = 80.41 GeV 
and oew = 1/128 |20|. And the other MSSM parameters are determined as follows: 



(1) For the parameters Mq jy ^ and Af h in sqark mass matrices, we set Mq = = 
= At = Ab = 200 GeV. 

(2) In the calculation of the corrections to the decays AP / XiXi : choose 
tan [3 and the masses of the two charginos m-+ as the input parameters to determine 

Al,2 

the relevant parameters in the chargino sector. Then the parameters ^ and M in the 
mass matrix of chargino can be extracted from these input parameters. (In this paper, we 
assume /x is negative.) 

(3) When we calculate the corrections to the decays A^ / X2X2' need four 
independent parameters to determine the neutralino sector. Usually these parameters are 
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taken as tan /3, M, M' and ^. But in our paper, we impose the relation between M and 
M': M' = (5/3) tan^ ^lyM |p and take tan/3, /i and m^o as the input parameters. 

In Fig. 2 and Fig. 3 we plot the relative corrections 6 as functions of tan /3 for ^4*^ and 
decays with niA = 250 GeV and niH = 250 GeV respectively. In the two figures, the solid 
curves represent the decays / —>■ xtxi with the input parameters = 110 GeV, 
m-+ = 300 GeV, and the dotted curves represent the decays /H^ X2X2 where we 
take m^o = 100 GeV, /x = —200 GeV. For the chargino decay of A^{H^) (solid curve), the 
relative correction 5 decrease from —7.5% (—13.2%) to —4.6% (—5.1%) as the increment 
of tan /3 from 2 to about 7, and then increase to —14.0% (—17.8%) as tan/3 to 36. The 
corrections to the chargino decays of ^ and are very sensitive to the value of tan /3 
and quite similar to each other. In contrast to the chargino decays of A^ and , the 
correction to the neutralino decay of ^ X2X2 (dotted curve in Fig. 2) is insensitive 
to the value of tan (3. It varies only from —5.4% to —7.3% as tan (3 running from 2 to 36. 
For the decay X2X2' there are two peaks at tan [5 ~ 1.85 and 30 which corresponding 

to the resonance effects at the positions where mH(250GeV)~ 2 x m^_^(2 x 124.65GeV) 
and m//(250GeV)~ 2 x m^^(2 x 124.07GeV) respectively. The absence of these resonance 
effects on the dotted curve in Fig.2 {A^ —>■ X2X2) due to the fact that the couplings V'^oftf . 
and V.Qrn, are zero when i = j (see Eq.( 5.20[ )). 

i J 

In Fig. 4 and Fig. 5 we present the relative corrections as functions of tjia for the chargino 
and neutralino decays of A^ respectively. Here tan /? is set to be 4 (solid curves) and 30 
(dotted curves) to make comparison. The values of the relevant input parameters are 
given in the figures. The peaks at niA ~ 348 GeV on the solid curves (tan /? = 4) in Fig. 4 
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and Fig.5 just reflect the resonance effect at tfie position ~ 2 x mj(2 x 174.3 GeV). 
As shown in the figures, the relative corrections are insensitive to niA and less than —7% 
for almost all the values in the range of 250(220) GeV to 400 GeV except for the region 
around the resonance position: ruA ^ 2 x mt- As mentioned in Appendix C, this top 
quark resonance effect are remarkable only when tan (3 is small for the appearance of cot (3 
in the coupling of the vertex — t — t. Therefore, there are no apparent resonance peaks 
at mA ~ 348 GeV on the dotted curves (tan /? = 30) . 

In the case of the chargino decay of — xtXi^ for tan/3 = 30, m-+ = 110 GeV 

and m-+ = 300 GeV (dotted curve in Fig.4), the squark masses are m^^ = 180.96 GeV, 
rrif^ = 322.28 GeV, mg^ = 241.57 GeV and mg^ = 160.71 GeV. The first squark resonance 
peak should appear at uia ~ nif^^ + m^^ = 402.28 GeV which is just beyond the upper 
limit of ruA (400 GeV), therefore it can not be displayed completely in Fig.4. As shown in 
Fig.4, the relative correction to this decay varies from —10.8% to —16.3% as increasing 
from 250 GeV to 400 GeV. 

For the neutralino decay of A'^: A'^ ^ X2X21 with the input parameters tan /? = 30, jU = 
—200 GeV and m^o = 100 GeV, the squark masses are m^^ = 179.65 GeV, m^^ = 323.01 
GeV, mg^ = 262.28 GeV and m^^ = 124.07 GeV. The sharp peak on the dotted curve in 
Fig.5 is just at niA = fn^^ +"^62 ~ 386.35 GeV. The correction increases slowly from —6% 
to —8% as the mass of AP varying from 220 GeV to 370 GeV, and then increases rapidly 
as niA to 386 GeV for the resonance effect around this position. 

Analogously, the mn dependence of the corrections to the chargino and neutralino 
decays of are depicted in Fig.6 and Fig.7. From Fig.6 (Fig.7) we can see that for the 
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chargino (neutralino) decay of there are three peaks on the curves. The peak on the 
solid curve corresponds to the resonance effect at the position tuh = 339.88 (315.14) GeV 
~ 2 X m^-^ = 2 X 170.27 (158.00) GeV in the case of tan^ = 4. The other two peaks on 
the dotted curve are just at thh = 321.06 (248.09) GeV ~ 2 x m^^ = 2 x 160.71 (124.07) 
GeV and mn = 361.05 (359.05) GeV ~ 2 x m^-^ = 2 x 180.96 (179.65) GeV. They can 
also be explained by the resonance effects of the chargino (neutralino) decay of H'^ with 
tan/3 = 30. Except for the resonance effects, the corrections are about —5% ^ —7% and 
insensitive to mn when tan/3 = 4. But, for tan/3 = 30, the corrections can exceed —10%, 
especially in the case of the chargino decay of H^, it can reach about —17%. 

The relative corrections to the chargino decays of A'^ and H*^ versus the lightest 
chargino mass m-+ with tan /3 = 4 and 30 are given in Fig.8. The upper two curves 
represent the cases of tan /? = 4 and the lower two, tan /3 = 30. The radiative corrections 
to these decay modes are all insensitive to m~+. In the case of tan/3 = 30, the relative 
corrections are considerable. They can reach about —11.7% and —15.2% for the chargino 
decays of and respectively. 

For the neutralino decays of A^ and H^, A^/H^ X2X2) P^o* relative correc- 
tions as functions of m^o in Fig.9. The curves in this figure are quite different in shape 
from those in Fig.8. The corrections are sensitive to m^o, especially in the case of large 
tan/3. For tan/3 = 30, there are two small peaks on the solid and dotted curves which 
are just at m^o = 128 GeV ~ m;, + m^^ = 4.3 GeV + 124.07 GeV. They are the only 
resonance peaks which can appear in this figure as m^o in the range of 100 to 160 GeV. 

In summary, we have computed the supersymmetric electroweak corrections to the 
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partial widths of the deacys XiXi and A^/H^ X2X2 in the MSSM. Only the 

third generation of quarks and squarks are included in our one-loop Feynman diagrams 
calculation. The corrections to these decay modes are all important and can exceed 
—10% in some parameter space, therefore they should be taken into account in the precise 
experiment analysis. 

After the submission of this manuscript to Physical Review D, we found similar calcula- 
tions in Ref.[23| where they gave the electroweak corrections to the decays {h^, H^, A^) ^ 



X° + X° and x° ^ (h^ H\ A^) + {h j = 1, 2, 3, 4). 
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5 Appendix 
5.1 Appendix A 

The matrices J^'^''^^^ and F^'^''-''^ are defined as follows: 

\ _ f cos a sin ct \ f Qij 
J^j J \ — sin a cos a J \ Sij 

-^ij \ = ( - cos/? sin/3 \ / Qij 

J Y sin/? cos/? y Y '^ij 

\ _ f — cos a sin Q \ / Q'-j 

f}j J ~ \ sina cosa J \ S^'j 

i^f A / cos/? sin/? \ / Q^;. 
F^j J \ - sin/? cos/? J V 



(5.1) 



16 



where 



Qij = 

Q'lj = 7;Wi3iNj2 - Nji tan9w) + Nj3{Ni2 - Na tanOw)] 



S'lj = \Wii{Nj2 - Nji tanOw) + Nj^{Ni2 - Na tanOw)] (hj = 1, ...,4) (5.2) 



2' 

Prom these definitions, we can find that Q", S" and F^'^'^'^ are symmetric matrices. The 
counterterms of the vertices A^{H^)xtxf ^^'^ A^{H'^)XiX'j can be written in the form 

+\H'fr{PLFlf* + PRFlf)x] - '{A'-X^iP,F^f^* - PnF^^)^ (5.3) 

where J^f^'^ and Pj^f^^'^ are the form factors of the counterterms. 
The expression of J^f^ is 



(5.4) 



^ k=i 

(c)A (c)H (c)A 

The other three form factors J-j ' , F^ ' and F^j ' can be written down analogously. To 
express the form factors F^^^'^^ and J^j^'^^ contributed by the vertex corrections, we 
introduce some notations to denote the vertices which will be used 

H^-t-t: Vhou H^-b-b: V^ou 

A'{H')-tl-ir. ^^o(^o),-|,-^. A'{H')-bl-br. ^^o(^o)StB, 

^-^-tr y^iPL + V^iPR b-^-b,: V^P, + V^Pn 

t-xt-h: V^iPL + %iPR b-^-i,: {V^,-P, + V^,^Pn)C{5.5) 

17 



For the decays AP /H^ — x^x°, we define 



Bq = Bo{q) = Bo{-p,mq,mq) 

Co%,i2 = ColiiM^J^Q^k) = Co,u,i2{ki,-p,mq^,mq,mq) 

Co?i,i2 = C'oSi,i2(^' J' 1^ ^' = Co,ii,i2(^i, -P, mq, mq^,mq^) (5.6) 



i^'^^^)=i^'^^^)(..,.,/c)^^o,,o,,l^,^V^o,, 
i^'-«=i^'-(^)(.,,.,.).F,o,,o,,^S,,^^^^ 

^ , V^^o,^o,,V^fo,,F4-, (5.7) 

where q = t or b, p = are the outgoing momenta of x° and x°. Then the 

form factors F"-^"^'^^ can be written as 



-iNr 



q=t,b 
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where Nr = 3 



k=l 
+ 
2 



"xr xr xr** ; ^ xy 

4^'^ = E h^^^'^^f ^ - -x?^/^^^^l? + (-,oF/(^) - r.,oFi^(^)) eg)] (5.10) 



k,l=l 



fc=l 

Af) = -..,oF5^)ci? + (n.,oF5^)-n.,o^ (5.11) 

fc,/=i 

For the decays /H^ —>■ xfxj^ we define 

Ml) _ Ml) (■ . „ = / C'o,ii,i2(A:*,-p,m^^,mt,mt), g = i 

"^0,11,12 — "-^0,11,121'' J' '''^ — 1 r> (h nTTim^n^ A 



(5.13) 



where /cj and kj are the outgoing momenta of xt xj- As in Eq.( p.7D and ( p.Sj) , the 
definitions of J^f^^^^ and J^^^^^"^^ can be obtained from the definitions of Ja'^^^^ and 
jA,H{2) ^^^y substituting for their superscripts L, L. Then the form factors J^j^'^^ can 
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be expressed as follows 



where 



J. 



v{A) 



v{H) 



(n 



H{1) 



+ ) \q=t ) \q=b 



(5.14) 



n 



^(1) 



Xi Xj ^ 



z 

fc=i 



n; 



A{2) 



kl=\ 



9^6 



+ m 



(5.15) 



n 



H{1) 



n 



H(2) 



z 

n 

k=l 



5n 



z 



jHi2)M2) 



Xi 



(5.16) 



2 2 

2 r 
mf =F , 



q = t 
q = b 



(5.17) 



In this paper, we adopt the definitions and calculation formulas of the two-point and 



three-point Passarino-Veltman integrals shown in Ref.|21] and |22|, respectively. 
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5.2 Appendix B 

The counterterm of the potential vjj^^^^ is 

+ ^H%\6ZHhmjj + 6ZhHml + 26bHh) + l^^^SbcG 

+ ^A''\6ZAAml + 6m\ + 6bAA) + lG^A\6ZAGml + 26bGA) (5.18) 

Here dbuHi^^hh-, ■■■ are functions of JT^/ and dT^. We can get their expressions from 



that of bHH,bhh, ■■■ in Eq.(2.16) by substituting 6Th and 6Th for Th- and Th- Then the 



renormahzed two-point Green functions can be obtained 

if^^^(p2) = i^p2_^j^)+ij:HH^p2^^-^^ZHH(.P^-ml)-6ml-6bHH] 

it^^ip^) = i{p^ -ml)+iT.>'\p') + t[5ZtaXp' -ml)-5ml-5by,h] 
it^\p') = iJ:^\p^)+i[^6ZHh{p'-ml) + ^6ZhH{p'-ml)-6bHh] 

if^^(p2) = i(^p2 _m\)+iJ:^^{p^)+i[6ZAA{p^ -m\)-6m\-6bAA] 

where the mass parameters of the CP— even neutral Higgs bosons mn and nih are deter- 
mined by mATrnz and tan/3 as in Eq.( p.l3| ). 
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5.3 Appendix C 

Here we list some of the couplings of Higgs-quark-quark and Higgs-squark-squark vertices 
which were used in the analysis of the numerical results: 

A -t-t. W =-^;;^7 

H'-i-t: F^o« = -zf^^ 

" " 2mwsml3 

^0 _ _ ^ . ^ ^^^.^^.^ =_£±.^i^ + A, tan /3) (7^;^7^52 - 7^f27^f 1) (5.20) 

(1) For large tan (3{> 30), the coupling Vao« is suppressed due to the factor cot /?, therefore 
the radiative corrections to the decays — > Xi^Xr/X2X2 contributed by the diagrams 
involving the vertex A^ — t — t are negligible. 

(2) The couplings V^oftj. and V^oft^. are zero when i = j. 
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Figure Captions 

Fig.l The relevant Feynman diagrams to the decays A^/H^ XiXp /H^ xfxj- 
(a) and (d) tree-level diagrams; (b) and (e) counterterms for vertices; (c-l)-(c-4) and (f- 
l)-(f-4) one-loop vertex corrections. In diagrams (a), (b) and (c), i = j = 2, (d), (e) and 
(f), i = j =1. In all the loop diagrams, the subscripts k and / run from 1 to 2. 

Fig. 2 The relative correction to A^ decay as a function of tan/3. The solid and dotted 
curves represent the decays A'^ —>■ xtxi and A'^ —^ X2X2 respectively. 
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Fig.3 The relative correction to decay as a function of tan f3. The solid and dotted 

curves represent the decays XiXi ~^ X2X2 respectively. 

Fig.4 The relative correction to the decay A° — xlxi as a function of niA- 
Fig.5 The relative correction to the decay X2X2 as a function of tua- 

Fig.6 The relative correction to the decay —>■ XiXi as a function of iuh- 
Fig.7 The relative correction to the decay H'^ X2X2 as a function of niH- 
Fig.8 The relative correction as a function of m -+ . The solid (dotted) curves represent 

the decay — xtxi ~^ xtxi) for tan/3 = 4 and 30. The input parameters are set 

to he mH = 350 GeV (via = 350 GeV) and m-+ = 300 GeV. 

Fig.9 The relative correction as a function of m^o. The solid (dotted) curves represent 

the decay — > X2X2 ~^ X2X2) tan /3 = 4 and 30. The input parameters are set 

to he mH = 350 GeV (m^ = 350 GeV) and n = -200 GeV. 
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Fig.l 
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Fig.7 
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